Route to nonlocality and observation of accessible solitons 



o 
o 

< 



C/2 

Q 

• I— I ' 
Oh 

q ; 
o 

^ ■ 



> ■ 

O 

oo . 
O ■ 

: 

o . 

^ : 



Claudio Conti, Marco Peccianti, Gaetano Assanto 
NOOEL - Nonlinear Optics and Optoelectronics Lab, 
National Institute for the Physics of Matter, INFM - Roma Tre, 
Via della Vasca Navale 84, OOI46 Rome - Ital-^ 
(Dated: February 2, 2008) 

We develop a general theory of spatial solitons in a liquid crystalline medium exhibiting a non- 
linearity with an arbitrary degree of effective nonlocality. The model accounts the observability of 
accessible solitons and establishes an important link with parametric solitons. 



In various areas of applied nonlinear science, non- 
locality plays a relevant role and radically affects the 
underlying physics. Some striking evidences are found 
in plasma physics^Q 0) 13 or in Bose-Einstein conden- 
sates (BEG), 0, I3 where, contrary to the prediction 
of purely local nonlinear models, nonlocality may give 
rise to, or prevent, the collapse of a (plasma- or matter- 
)wave. In nonlinear-optics, particularly when dealing 
with self-localization and solitary waves, nonlocality is 
often associated to time-domain phenomena through a 
retarded response (see e.g. spatially nonlocal ef- 

fects have been associated to photorefractive Ifll Hfl 
and thermal or diffusive responses To assess 

the role of nonlocality, theoretical studies tend to dis- 
tinguish between highly and weakly nonlocal behaviors 
by comparing the spatial extent of the ma- 
terial response (the so-called kernel-function) and the 
optical beam waist. Specific kernel functions, however, 
strongly depend on the physical system and, as in the 
case of BEG, .31 they are hard to determine and apply 
to experimental results. On the other hand, they are at 
the basis of the theory of spatial optical solitons (SOS) 
in highly nonlocal media. 

SOS have become the subject of intense theoretical and 
experimental investigations, both on the grounds of their 
packet nature and in view of applications, particul arly 
in the exploitation of their wave-guiding character. |l7l | 
Diverse material properties have been studied in conjunc- 
tion with SOS existence and properties, including various 
mechanisms able to counteract diffraction in one or both 
transverse dimensions. 0, Spatial solitons due to 
a local nonlinearity have been known since the original 
work of Ghiao, Garmire and Townes with reference to 
Kerr media. 20] In 1997 Snyder and Mitchell 1^] inves- 
tigated SOS in a highly nonlocal system, i. e., a medium 
exhibiting a power - rather than intensity - dependent re- 
fractive index. They introduced the term accessible soli- 
tons for those spatial solitary waves, owing to the sim- 
plicity of the theory and transverse profiles obeying the 
2-dimensional equation of a quantum harmonic oscillator 
(which gives a gaussian profile like the so-called "gaus- 
sons" 21]). Shen pointed out to this connection as an 
intriguing one between distinct fields of modern physics 
and, provided the required large correlation lengths could 



be made available, to the demonstration of accessible soli- 
tons as a challenge well worth undertaking, [iij 
In this Letter we introduce a model able to describe op- 
tical spatial solitons and the smooth transition from the 
purely local (in the limit of a Kerr nonlinearity) to the 
entirely nonlocal case. While pursuing a general the- 
ory, however, we chose to address a specific and available 
nonlocal system, i. e., nematic liquid crystals (NLG) in 
a planar cell. NLG have been proven to exhibit a sub- 
stantial nonlocal nonlinearity of molecular origin .23.1 and 
to support (2-fl)-dimensional spatial solitons, l24l even 
in the case of spatially incoherent excitations. 25] Af- 
ter deriving the ruling equations and defining a suitable 
nonlocal parameter able to span the soliton-family from 
pure-Kerr or Townes-like (T-) to higly nonlocal or acces- 
sible (A-) solitons, we will outline the rather remarkable 
connection between our model and the equations describ- 
ing quadratic two-color solitons (or simultons) in para- 
metric media. [2^ 

Let us consider the simple geometry sketched in figure 
1 : a planar glass cell containing an undoped NLG with a 
preset orientation of its molecular director. The aligned 
liquid crystal, anchored at the bounding interfaces, be- 
haves as a positive uniaxial with n//=nz and nj_=n^. 
In the presence of an external quasi-static (electric or 
magnetic) field or special anchoring at the interfaces, the 
refractive index n{0) in the (x, z) principal plane can ex- 
hibit a distribution along x, with n± < n{9{x)) < n//. 
|27j For a light beam linearly polarized along x and prop- 
agating along z, with transverse size well below the cell 
thickness L, neglecting vectorial effects and adopting the 
paraxial approximation, the evolution of the optical en- 
velope A is described by the Foch-Leontovich equation: 
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^\yA kln\\sin{Bf - sm(6'o)^]A = (1) 



where feg is the vacuum wavenumber. 
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optical anisotropy, fc^ = fco(n^ -I- n^sm(0o)^), Q the tilt 
angle of the NLG director, and the tilt in the absence 
of a light beam. When an external electric field applied 
along X and an optical excitation as in equation (1) are 
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present, 9 is subject to re-orientation according to [2 



applied voltage 
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sin{29)+ " °' ' sin{29) = 

(2) 

with K the relevant elastic constant taken equal for splay, 
bend and twist, Aenp the low-frequency anisotropy, and 
E the rms value of the quasi-static field. 
In the absence of a light wave, therefore, the orientation 
angle 9 is determined exclusively by E and, due to sim- 
metry, only depends on X:[23j 
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sin{2d) = 0. 



(3) 



In our case the boundary conditions correspond to the 
planar alignment: 9{X = -i/2) = 9{X = L/2) = 0. In 
the general case, the angle distribution can be written as: 

9{X,Y,Z)^9{X) + ^^{X,Y,Z). (4) 

Taking the cell much larger than the beam waist, we can 
use and Q in Q and neglect the derivative of 9. For 
9 = 9f), at the first order in 5* we obtain: 



K^'XYZ 



XY 

^1 — i^iRpE^-q, 



(5) 



having chosen 0n — 7r/4 in order to maximize the nonlin- 
ear response. |27| We write Q in a dimensionless form 
by setting A = {Ac/a)a{R/Rcy/a, Z/aZc)exp{iZ/aZc), 

= {'^c/a)i:{R/Rcv^,Z/aZr) with A^ 
SAejipE^/n^oklnlK, Z, = 2kRl,Rl = TrK/2AeRFE^, 
^-c = 2AeRFE'^/TTklnlK, and e = AtRpE^ /2TTk'^K . 
a is a free parameter, to be used hereafter to 
trace the family of spatial solitary waves, and 
{x,y,z) = [Xj ^/olRc^Y I \foLRc^Z I aZ^ are normal- 
ized coordinates. The resulting system is: 



; d a 



(6) 



To underline the physical meaning of the parameter a, 
let us consider the case e = and formally write the 
solution of the re-orientation equation as ■0 = (1 ~ 
V^/a)^^|ap/(2a). For large a, we obtain for the field 
envelope: 
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Equation lO rules collapse-free weakly-nonlocal 
media, 3], while the neglecting of terms such as 0(1/q;^) 
describes light propagation in Kerr media which, on 
the contrary, are subject to catastrophic self-focusing. 
Reducing a increases the degree of nonlocality in the 
interaction between medium and optical beam. In the 




FIG. 1: Planarly aligned nematic liquid crystal cell for the 
observation of spatial solitons. 9iX) is the low-frequency 
voltage-induced angle displacement. 8o = 9(0) and ^ is the 
perturbation due to the propagating optical field. 



following wc will employ a as an arbitrary parameter 
spanning the whole family of SOS; we expect large a 
values to be associated to T-solitons, whereas small a 
will address A-solitons. Note that, according to 
when e = the fundamental solitary wave solutions 
of Eqs.© are stable because they realize an absolute 
minimum of the Hamiltonian. 

Solitary solutions of 10) are defined by dz = 0. With- 
out loss of generality, taking a real-valued we have: 



V^a -a + atp = 



(8) 



Noteworthy enough, system © is identical to what de- 
termines the profile of parametric spatial solitary waves 
in x^^^-media. [ifllE^ This makes an unexpected connec- 
tion between self-trapped beams in two distinct physical 
systems encompassing rather diverse nonlinear optical re- 
sponses: the ultrafast electronic nonlinearity of quadratic 
crystals and the slow molecular re-orientation of liquid 
crystals. The similarity is only limited to the profiles of 
the solitary waves, while their dynamic properties, such 
as stability, are quite different. As shown below in a 
practical case, to a first approximation e is negligible, 
thus the angle "adiabatically-foUows" the light beam; in 
general this is not true for the harmonic-field of a para- 
metric soliton. It is well known, in fact, that for small 
a the V'-field (the second-harmonic for x^^^ crystals, and 
the re-orientation for NLC) is much wider than the a- 
field (see Torruellas et al. in The opposite holds 

true for a — s- oo: in the x^^''"SOS literature, this is the 
Kerr-limit, its dynamics resembling that of x^^^ materi- 
als. In NLC, for optical beams much wider than the re- 
orientation profile, T-solitons approximate well the solu- 
tion (in the framework of the validity of the large cell ap- 
proximation) . The angle-profile perturbation is localized 
close to the optical beam axis, as typically pointed out 
when addressing the Kerr-hke response of NLC. [i^l 
Here we are rather interested in the opposite limit. 
In figure 121 we graph the numerical solutions of eqs. © 
as obtained by a relaxation procedure. The ratio p be- 
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FIG. 2: Ratio p between the waists of the angle tp and the 
field a from eq.l|HJ. The inserts show two different profiles 
(dashed line: solid line: a) for values of a addressing the 
two hmits in the family of solitary waves. 



tween SOS beam- and angle- waists (standard deviation) 
is plotted versus a. When a approaches zero, i. e., for a 
response-length extending well beyond the optical waist, 
the V-'-field is much wider than the a-field. This is the 
high-nonlocality regime. In the figure we indicate the 
two opposite limits, namely A- and T-solitons. In the A- 
limit, however, it is worth proving that eqs. © reduce to 
the model in [lj| . The authors in 14] postulated a highly 
nonlocal medium in which the index of refraction could 
be expressed as = ng — a^{'P)B?^ with V the optical 
power. Here we show that the NLC re-orientational non- 
linearity does indeed exhibit such feature, as speculated 
rather skeptically in ji^ . 

To study the highly non-local regime (a 0), first 
we solve eqs. (jSj), in the regions around r — Q {r = 
\/ + 2/^) and r — s- oo, and then match the resulting 
expressions. For ?- = we introduce the expanded vari- 
ables 113 C = r/a^/^ and ^ = r/a^/"^. Next, we ex- 
press the field as a = --^/(^) with /(O) = 1, while 

^ = V(°HC) + (l/Va)V'^"^nC)- At order 0{l/a^/^) we 



obtain for ij): 



V'J? + {l/O^T + (l/2)aV'(«'^\) = (9) 



For a ^ with Q fixed, the solution reads i}}'"^^ = ipo — 
a = V'o - a^C^/Sy^. At order 0(1/q;^/^) we have 
= 0; at order 0(l/a^/^), the field equation reduces 



^2 2 

a r 



to the harmonic oscillator equation in /: 

hi + + ^^"'V - ima^f = (10) 



corresponding to the result in |l4'|, and V'o = 1 at or- 
der 0(1). After some algebra, the perturbative approach 



provides the SOS profile near the origin (?/;^ 



1) 



V' = 1 + 



J, 2 



2/^0^): 



(11) 



with ^0 an arbitrary parameter. The appearance of an- 
other free parameter refiects that, when a = 0, eqs. © 
are invariant with respect to the transformation {x,y) — )■ 
{x / / ^),a /ia,V' ^ I — ^ + with fj, arbitrary. 
The physical meaning of Hll|) can be elucidated by ex- 
pressing them in terms of the normalized beam power 
P = J J a^dxdy. a = (%/2P/47r)ea;p(-r2p/167r), and 
■0 1 + (P/47r) - (P/87r)r2. The latter results can 
be also obtained by using a multiple scales expansion 
in -y/Pr, Pr, .... For the field intensity in the original 
variables(P = Va^+T^), we have: 
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■exp[ 



(12) 



In H12() Zq and rj = Z^jn are vacuum and medium 
impedances, respectively, c the speed of light in vacuum, 
n = fc/fco = \J n± + 71^/2 when ^0 — '''/4, and Vc a ref- 
erence power which depends on material properties and 
cell-polarization: Vc = IGcnAeRpE^ / kln^. The soliton 
profile is gaussian and completely determined by V. Ren- 
dering explicit the relation between V and the (intensity) 
waist W, its existence curve is: 



W2' 



(13) 



Let us now investigate the region of large r. Accord- 
ing to eqs. © the asymptotic behavior of t/j as r — s- 00 is 
given by the modified Bessel function: t/^ — > CKQ^-Zar), 
governing the angle decay far from the beam axis (G is 
a constant to be determined). The angle in Hll() and its 
derivative must match the expression at infinity, provid- 
ing the turning point tt between the two regions: in the 
presence (r = 0) and in the absence (r 00) of the 
optical excitation, respectively. We end up with 
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4^ '"^^^8^' 



P ^ 2rTKo{^r) 



(14) 



which can be further simplified by taking into account 
that, for large arguments, the Bessel functions Kq and 
Ki have the same asymptotic behavior. The approxi- 
mated solution is = \/a{8n / P); being TZt — y/aRcrr 
it reads: Ut /TZc = Vc/V. When V » Vc the profile of 
the angle is dominated by the modified Bessel function, 
and decays with typical length Re- Therefore, it seems 
natural to label "highly nonlocal" a regime when the an- 
gle profile is much larger than the beam waist, such that 
TZc » W- When the power is much greater than Vc 
or, equivalently, when the soliton waist is much smaller 
that the extent TZc of the elastic response, we are in the 



A-soliton regime. The power dependent perturbation of 
the refractive index is: 



(15) 



with ric — AeRpE^ /T^nkQK . Moreover, our approach 
enables to go beyond the harmonic osciUator. By in- 
troducing the new transverse scale a — rjc?!"^ we can 
solve the resulting equation for ■(/''°-' at order 0(a^l'^\. 

+ (l/cr)V'^''^ - 4ctV^§ = 0. After some algebra 
we obtain the new approximation in normalized units: 
Vj = 1 + (P/47r) - (P/87r)r2 + [P j^-nfr'^ j \. Then, the 
corresponding power-dependent refractive index pertur- 
bation reads: 



(16) 



Thus higher-order approximations imply anarmonicity of 
the nonlocal potential, with a refractive index still de- 
pending on power. As the latter increases, higher powers 
of the ratio V jVc must be taken into account, similar to 
local media with regards to the powers of the intensity. 

Finally, comparing our theory with an actual experi- 
mental geometry, such as employed in [i^], typical pa- 
rameters for a 514nm wavelength and the E7 NLC (in SI 
units) are: ^ 1, K ^ 10"", £; = 1.3 x 10-'^,L = 75 x 
10"^ AcRF = 20 eo. Correspondingly, "Pc = 2 x 10~'^W, 
Rc = 22/im, and e — 10"^. Since V inside the cell was 
of the order of O.lmW, we may state that in experi- 
ments the highly nonlocal regime {V >> Vc) is being 
addressed. The soliton waist, after is of the order 

of 3/iTO and in agreement with the reported results. Note 
that An(0,P) = 5 x 10~^, and the angle perturbation ^ 
of the order of 10~^ radians, thus justifying the adopted 
model. 

In conclusions, for the first time to our knowledge, 
we have presented a self-consistent analytical theory of 
two-dimensional spatial solitary waves in nonlocal me- 
dia. Our model has an intriguing unifying character, as 
it embraces several physical systems in which light self- 
trapping has recently been investigated. We believe that 
most of the observed SOS' in nematic liquid crystals are 
indeed accessible solitons, in-as-much-as NLC are highly 
nonlocal. This shines new light on self-localization in liq- 
uid crystals. Furthermore, we have presented the first 
derivation of a power-dependent constitutive relation for 
a real physical system, never reported elsewhere. We con- 
fide that our results will stimulate new experiments to- 
wards a deeper understanding of self-trapping in (highly) 
nonlocal nonlinear media and the development of novel 
all-optical devices, fioj 
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